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Abstract. The contribution to the entanglement of an impurity attached to one end 
I' , of a J — J2 ^ <^ quantum spin chain {S = 1/2) is studied. Two different measures of 

?-H ' the impurity contribution to the entanglement have been proposed: the impurity- 

^ ■ entanglement-entropy S'imp and the negativity M. The first, S'lmp, is based on a 

subtractive procedure where the entanglement-entropy in the absence of the impurity 

is subtracted from results with the impurity present. The other. A/", is the negativity of 

a part of the system separated from the impurity and the impurity itself. In this paper 

' '^ \ we compare the two measures and discuss similarities and differences between them. 

In the J — J2 ~ S model it is possible to perform very precise variational calculations 
Q I close to the Majumdar-Ghosh-point ( J2 = •//2 and (5 = 0) where the system is gapped 

with a two-fold degenerate dimerized ground-state. We describe in detail how such 
PsJ ■ calculations are done and how they can be used to calculate Af as well as S'lmp for 

^ I any impurity-coupling Jk- We then study the complete cross-over in the impurity 

entanglement as Jk is varied between and 1 close to the Majumdar-Ghosh-point. In 
^^ , particular we study the impurity entanglement when a staggered nearest-neighbour- 

C\| ' interaction proportional to 6 is introduced. In this case, the two-fold degeneracy of the 

f~^ ■ ground-state is lifted leading to a very rapid reduction in the impurity entanglement 

^3 ' as 5 is increased. 
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1. Introduction 

Entanglement in quantum spin cliains lias received considerable attention. For gapless 
quantum spin chains detailed predictions [1] of many aspects of the entanglement have 
been obtained from conformal field theory (CFT). For these (1+1) dimensional models a 
precise understanding of the entanglement has proven invaluable for the understanding 
of numerical techniques such as as the density- matrix renormalization group |2] (DMRG) 
as well as for the development of new techniques. The contribution to the entanglement 
arising from impurities has also attracted considerable interest [H IU |5l O [71 El ISl [TOl 
ITTl [121 [13]. Here the term impurities is used in the general sense and includes the effects 
of boundary magnetic fields [3|, boundaries |4J, qubits interacting with a decohering 
environment O [6l [71 [8] as well as Kondo-like impurities [9l [TOl [HI [121 [13] • The presence 
of the impurity can lead to different conformally invariant boundary conditions which 
can dramatically alter the entanglement. For a review see reference [H]. Quantum 
spin models have been used for studies of qubit teleportation and quantum state 
transfer [151 (HI (13 [HI [El [201 1211 [221 [23], where the change in entanglement arising 
from the impurities plays a crucial role. 

Usually the entanglement is defined in terms of the von Neumann entanglement- 
entropy of a sub-system A of size / and reduced density-matrix p^ defined by jMl [25] , 

5(/,L) = - Tr [pA In pa] , (1) 

where L stands for the total system size. It is known [4] that the entanglement-entropy 
of Heisenberg-spin-chains with open boundary conditions have a uniform as well as an 
alternating part: 

Sil,L) = Suil,L) + i-iySa{l,L), (2) 

where the subscripts u and a stand for uniform and alternating respectively. 

Based on this observation, the impurity entanglement-entropy Sinip{JKj,L) was 
introduced to quantify the contribution an impurity has to the entanglement of a spin 
chain j9l ^U\. To model an impurity, the interactions on one site are scaled by a factor 
of Jk- The impurity entanglement-entropy is given by the difference of the uniform 
part of the entanglement-entropies of the chain with the impurity present and the chain 



^ / + 1 L 
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Figure 1. An example of a spin chain with only nearest- neighbour-interactions and 
an impurity on the first site. The interactions of the impurity site are multiplied by a 
Jk- The chain is separated into part A (grey box), which includes the impurity and 
has the length I, and the rest of the chain. 



Impurity Entanglement in the J — J2 — 6 Quantum Spin Chain 3 

without the impurity. The sub-system A in the definition of the entanglement-entropy 
in equation [1] is chosen to contain the first / sites of the chain and thus contains the 
impurity site, too (see figured]). Upon removing the impurity site, A as well as the whole 
chain are one element shorter. The uniform entanglement-entropy in a system where 
the impurity is absent is thus given by Su{^,l — 1,L — 1). This leads to the following 
definition of Si^ap'- 

SimpiJx, ^, L) = ^^(with impurity) — ^^(no impurity) 

= 5„(J^,/,L)-5„(1,/-1,L-1), l>l. (3) 

This definition is similar in spirit to experimental procedures for extracting impurity 
contributions to, for instance, susceptibilities. It is also possible to consider the 
alternating part of the impurity entanglement but we shall not do that here. 

However, it was later pointed out [TT| [T2] that a more consistent definition of the 
impurity contribution to the entanglement can be obtained from the negativity |26j . 
The negativity between a sub-system A and the rest of the system is defined by 

A/-=^^. (4) 

where the Aj are the eigenvalues of the partial transpose of the density-matrix with 
respect to either A or the rest of the system. In contrast to the entanglement-entropy 
it can also be used to quantify the entanglement if the system is not in a pure state. 
This makes it possible to use the negativity to directly quantify the entanglement of the 
impurity-spin and another part of the system. 

It is an important question to what extent these two quantities yield the same 
information about the impurity entanglement. In this paper we compare these two 
measures and show that they agree for many universal features but are fundamentally 
different with M being the more general measure of the impurity entanglement. 

Entanglement and in particular impurity entanglement in gapped quantum spin 
chains has received comparatively little attention. In part this is due to the fact that for 
non-critical spin chains with a finite correlation length S, one expects [T] for a partition 
into 2 semi-infinite chains: 

5* ~ log ^ . (5) 

For instance, for the S = 1 Affleck-Kennedy-Lieb-Tasaki (AKLT) chain [271 [28] where 
the entanglement can be calculated analytically |29| |30| [3T| [32| |33| 131] 135] it is known 
that S{1, L) approaches a constant for large /, L in an exponential manner [Ml l35] on a 
length scale equal to the bulk correlation length ,^ = 1/ ln(3). At the Majumdar-Ghosh- 
point (MG-point) [36] of the J — J2 spin chain, where J2 = J/2, one finds that S{1, L) is 
either or ln(2) for any /, L. However, the impurity contribution to the entanglement 
can still he long-range in gapped spin chains. For the J — J2 spin chain this was shown 
to be the case at the MG-point [9l ^U\. This is due to the two- fold degeneracy of the 
singlet ground-state, corresponding to the two different dimerization patterns, and the 
associated presence of solitons separating regions with different dimerization. The model 
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we study is the slightly more general J — J2 — S quantum spin chain defined as: 
H = Jk[{1 + 5)So-Si + J2So-S2 

L-3 

+ 5Z [(1 + {~sy)s^ ■ s^+i + J2S,- s,^ 



i=l 



+ (1 + {-6)'^-')Sl^2 ■ Sl^i , (6) 

where we implicitly have set the nearest-neighbour-coupling J = 1. The coefficient Jk 



Figure 2. The J — J2 — 6 chain with an impurity coupled with strength Jk- 

describes the coupling of the impurity-spin and 5 is a staggering of the nearest-neighbour- 
coupling inducing dimerization. We focus exclusively on chains with open boundary 
conditions and the impurity at one end of the chain. With 6 = this model undergoes 
a transition from the gapless Heisenberg- phase to a dimerized state at J| = 0.241167 
|37ll38] . In the Heisenberg-phase a correspondence to the low-energy physics of a Kondo 
impurity can be established [39l HO] . This correspondence becomes exact at J|. The 
mapping to the Kondo-problem requires one to use open boundary conditions and we 
therefore do not consider periodic boundary conditions here. However, entanglement in 
models similar to the one considered here but with periodic boundary conditions have 
also been studied, see section 4 of Ref. jl4] for a review. 

With periodic boundary conditions and an even number of sites, the system at 
the MG-point has an exactly known two-fold degenerate dimerized ground-state and a 
gapped spectrum [361 l¥T] . (With open boundary conditions and L even the ground- 
state is unique with the same energy). The lowest lying excitations of the MG-chain 
are pairs of unbound solitons and upon introducing a finite dimerization the solitons 
become confined |l2]. For a chain with an odd number of sites solitons separate regions 
with the ground-state dimerization-pattern of the chain with an even number of sites 
[in 142] . In the gapped phase, in the proximity of the MG-point, it is possible to perform 
very precise variational calculations ^T\ H3| HI] for the J — J2 model (corresponding to 
(5 = in this study) since appropriate states spanning the variational sub-space can be 
identified with relative ease, as the ground-state at the MG-point is known. Within this 
variational framework it is possible to evaluate the entanglement ^U\. Here, we explain 
in detail how similar calculations can be performed for the negativity and discuss the 
appropriate sub-spaces used in these computations. Along the so called disorder-line, 
where J2 = (1 — S)/2, one of the two degenerate ground-states of the Majumdar- Ghosh- 
chain remains the exact ground-state. Close to the disorder-line, one therefore expects 
the same variational methods to yield very precise results. 
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In the past, these variational calculations have been performed more or less by 
hand [1T| H3| |U]. If the variational problem is, however, formulated as a generalized 
eigenvalue-problem, it is not difficult to automate the calculation, as for any two states, 
(fii and ipj, with all but no, one or two sites being in diniers, {(pi\SiSj\(pj) as well as the 
overlaps between different ip can be computed in an automated fashion. One advantage 
of the automated computation of these matrix-elements is that one avoids lengthy and 
error-prone calculations. Secondly, solving a generalized eigenvalue-problem can be 
done using standard libraries. This approach is not limited to the specific Hamiltonian 
studied in this paper. It can in fact be used to perform variational calculations using 
dimerized states forming the subspace in which the diagonalization is carried out for 
any Heisenberg-spin-Hamiltonian. 

In addition to providing a more direct analytical insight into the physics, the 
main advantages the variational approach has over other methods such as the DMRG- 
technique are computational cost and simplicity. This advantage does, however, not 
come for free. To identify a good subspace to use for variational calculations, one 
always has to rely on knowledge gained obtained by other means. This role is often 
played by the DMRG-technique. 

The paper is organized as follows: In sections [21 [3l [Hand |5] the variational approach 
is discussed along with the different sub-spaces used in the calculations. Section [6] 
describes the calculation of the negativity once the variational ground-state is known 
and sample results for the negativity are compared to DMRG-data. In section [7] the 
negativity and the impurity-entanglement-entropy are compared. Section [8] presents 
our results for the negativity and 5'imp for general Jk at the MG point. In section [9] 
we describe how the impurity entanglement is affected by the presence of a non-zero 
explicit dimerization, 6. 

2. The Variational Calculation as a Generalized Eigenvalue-Problem 

For completeness we review here how the variational problem can be formulated as 
a generalized eigenvalue-problem. Consider a set of states {</?«} to be used for the 
variational calculation. We then wish to minimize 

where the trial-state ip is given by: 

The minimization is to be performed with respect to the Cj. In variational calculations 
we make use of the fact that: 

'^^■- Eo< ^^^ , (9) 
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where Eq is the ground-state energy. One strives to find a well chosen set of (fiS as well 
as the state ipmin ^ span{y9j} that satisfies 



{(p\(p) v™„ (v^lv^) 

This state is the variational estimate for the ground-state. Let us denote the energy- 
eigenbasis with {Cj}. The variational state (f can evidently be written as (f = "^aiti. 
Thus 



E' 



EJ«iP 

> E, , (11) 

where E^ is the energy of the eigenstate c^ of non-zero projection onto span {ifi} that 
has the lowest energy. With any orthonormal basis {b,} of spanjy^j}, the variational 
guess for the ground-state ipmin can thus be found by diagonalizing H: 

• Minimization problem = Finding lowest eigenvalue of {bi\Hbj) 

If we do not have an orthonormal basis at our disposal, we have to do a little 
more work. Let us suppose that the ifi are linearly independent and the Cj are the 
orthonormal basis of energy-eigenstates. For the variational ground-state given by 
'^min = Yli ^iVi ^"^^ the Variational estimate of the ground-state-energy Ey we know 
from the equations (fTTj) that 

H N^ ajLfj = Ey 2, o,j^j ■ (12) 

j j 

We will now show that a vector a = (ai, . . . , a„) is a solution to the variational problem 
if and only if 

Wi : y^^aj{ipi\H(fj) = Ey^aj{ipi\ipj) . (13) 

j j 

With the definitions 

S}ij = {ipi\H(fj) and ^ij = {fi\fj) , (14) 

one can write equation (IT3l) as 



Sja = Ey^a, (15) 



which defines a generalized eigenvalue-problem. That equation ( TT2l) implies 
equation ( IT3|) can be seen by taking the scalar-product of both sides in equation (IT2|) 
with (fi. 

That equation (IT3l) implies (IT2l) will now be proved. We do this by proving that 
for any basis (linearly independent set of states that span the space) {(fi}: 

\/i : (ipilq) = {ipi\x) =^ q = t , (16) 
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for any vectors q, r. With q = HJ2j <^jfj and x = E^ ^ . aj(fj the relation flT6ll becomes 
the relation equation (fT3|) =^ f fT2|) . 
Proof: 

V^: (^,|q) = (^,|r) (17) 

V^: ((^i|q-r) = (18) 

=^ 

Wz: ^i±{q-t) (19) 

=^ 

q - r = D (20) 

We have thus proven that the minimization-problem of equation (ITO|) can be 
replaced by the generalized eigenvalue-problem of equation (fT5|) . if the states used to 
span the subspace in which the minimization takes place are linearly independent. 

3. The Hilbert Spaces Used in the Calculation 

The implementation of the variational technique as a generalized eigenvalue-problem 
allows us to just solve a linear algebra problem instead of solving a set of non- linear 
equations which is necessary if one approaches the minimization naively. The one 
requirement for this implementation to work is for the states that span the Hilbert- 
space in which the minimization is done to be linearly independent. This flexibility 
makes it possible to use different bases according to the speciflc needs of the problem. 
Here we are concerned with the two cases of even and odd length chains. 

(i) Chain with odd number of sites: 

The simplest ensemble of variational states for a chain with an odd number of sites 
are the so called single soliton states |4H H3| HU [10]. These states can be generated 
by taking L — 1 of the sites of the chain to be maximally dimerized while leaving the 
left over site (the soliton) to be in the 5*^ = 1/2 state. Here and in the remainder 
of the paper we use L to denote the length of the chain. The number of dimers 
between the undimerized site and the left end of the chain can be used to label the 
states. In flgure|3]we show a pictorial representation of the flrst three states. When 
drawing a dinier we use an arrow to denote the order in (| fj,) — | l.'\))/\/2. For an 
odd length chain with L spins there exist A^ = (L + l)/2 such states. Note that we 
here do not include states with the soliton in the Sz = —1/2 state. 

(ii) Chain with even number of sites: 

In studies of the entanglement of the impurity-site with the rest of a chain with 
an even number of sites, a suitable set of states can be deflned by leaving the bulk 
[L — 2) sites of the chain dimerized while one of the spins is in a valence bond 
with the impurity. Example-states are given in flgure HI These states are chosen 
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♦ — — -^ 

(a) (fo 

•— •— • f •— >— • • > • 

(b)<^i 

•— — • •— — • f •— •— • 

(c) ^2 

Figure 3. The first three variational states used in the calculations with the single 
soliton space for a chain with an odd number of sites. Note the arrows on the dimers 
that indicate their direction. 




(c) </?2 



Figure 4. The first three variational states used in the calculation of the negativity 
for a chain with an even number of sites. 



in order to reflect the presence of an impurity valence bond (IVB) [10], connecting 
tlie impurity site to tlie bulk of the chain. 

4. The Matrix Elements of of the Hamiltonian S) 

To set up the generalized eigenvalue-problem it is necessary to calculate the matrices 
S^ and 03 which were defined in equation (IT^ . In this section we will explain how to 
compute Sj. In order to allow for the most general Hamiltonian possible we write it as 
follows: 

H = 2_^ JmlSmSl ■ (21) 

{m,l} 

Here, the summation {m, 1} is over 'bonds' of the Hamiltonian and the coupling 
constants J^i can take any value. It is convenient to write the matrix-elements S)ij 
in terms of the operators h^i = SmSi — 1/4, because they have a simple action on the 
(fs. This gives: 

{(Pi\H(fj) = ^ Jml{fi\SmSllfj) 
{m,l} 

= 5^ Jm/(</'i|[^m« + l/4]^j) , (22) 

{m,l} 



Impurity Entanglement in the J — J2 — 6 Quantum Spin Chain 



The following description of this action is an extension of results presented earlier 
in the context of quantum Monte Carlo simulations in the valence bond basis where it 
is necessary to take the hmi to have the opposite sign of what we use here. The action 
of hjni on states relevant to us is given by: 



hmi[m; l][k; n] 
hmk[m;l][k;n] 

hmk[m;l] tfc 



— [m; l][k; n] 
-[m;k][n;l] 



tm [/; k] + [m; I] Ia 



(23) 
(24) 

(25) 

(26) 



where 



h;^] := -^( tmifc - t4m ) 

are the spins at m and A; in a singlet-state. 

In figure \5\ a pictorial representation of the equations 023p - (l26p is shown. To fix 
the phase of the states we include in the calculation, we again represent [m; k] by an 
arrow pointing from m to k. For the calculations we present here these are the only 
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Figure 5. The action of hmi on the variational states. The dotted line signifies the 
sites acted upon by hmi- Arrows fix the phase of singlets between sites of the lattice. 
The convention is explained in the main text. 



rules needed. However, when considering the action of the Hamiltonian it is sometimes 
convenient to apply the simple rules 

[m;k] =-[k;m] (27) 

[k; I] U = [k] m] ti + tfc [m; I] , (28) 
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in order to reduce the action of the h^i to one of the above. 

We want to stress that even though the single sohton states can be orthogonahzed 
with relative ease jl6] this is not generally the case and even after orthogonalization the 
resulting eigenvalue problem is still non-trivial for a general set of Jim- Whether or not 
it is possible to proceed with analytical calculations often depends on specific choices for 
the coupling constants in the Haniiltonian. Here, we can solve the variational problem 
for any values of the couplings. 

As can be seen in equation ( l24l) and figure [5]^b), the application of hmi to a given 
state ipi may result in a state that is not element of span{(^j}. Therefore, in order to 
evaluate 

{(Pi\Hipj) = ^ Jmli'^ilhml'^j) 
{m,l} 

we need to compute the overlap of hmifj with all the variational states. This slows 
down our calculations since it is not sufficient to build a table of all the overlaps of the 
(fii. We now turn to a discussion of how the overlaps are calculated. 

5. The Matrix Elements of the Overlap-Matrix ® 

As is explained in [171 H5] the magnitude of the overlaps of valence-bond-states depends 
solely on the length of loops created by overlaying the two states. If one chooses I to be 
the length of such a loop then the overlap is of magnitude 2^'/^+^. To get the sign of 
the overlap we, in addition to overlaying them, reverse all arrows on one of the states. 
Then we follow the loop and count how often one has to go against the arrow. If the 
resulting number is odd, the sign of the overlap is negative. If the resulting number is 
even, the sign of the overlap is positive. Because the total number of bonds in the loop 
is even this leads to a well defined sign of the overlap. To get the overlaps of two states 

(a) 




(c) 

Figure 6. The two states (a) and (b) together with tlie corresponding loop-structure 
created by reversing the arrows of one state and overlaying the two states (c). Both 
loops contribute a minus-sign. The overlap is given by 1/4. 

one only has to follow this prescription for every loop and multiply contributions from 
different loops. Figure M^c) shows the loop-structure for the states shown in figure M^a) 
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and figure El^b). The left loop is six bonds long and going around it one has to go against 
the direction of the arrow three times. The right loop is 2 bonds long and one has to 
go against the direction of the arrow once. The overlap is therefore given by: 

(-2-^)(-2°) = \ . 

In the case of unpaired spins being part of the states the rules have to be modified 
slightly. If there is one unpaired spin in both states, upon overlaying the two states, 
there will in addition to loops be a string that goes from one unpaired spin to the 
other. As is the case for the loops, for the string the magnitude of the contribution to 
the overlap depends only on the length. One finds the contribution to be: 2~'/^. The 
contribution to the sign of the overlap can be determined in exactly the same way as in 
the case of complete loops. 

6. The Calculation of the Reduced Density-Matrix and the Negativity 

For the calculation of the negativity describing the entanglement between the impurity- 
spin and a distant part of the chain, it is necessary to separate the chain into three 
parts [26] [TT] : the impurity-spin, the part whose entanglement with the impurity-spin 
is to be quantified and the part that separates the two (the parts will in the remainder 
of the paper also be referred to as region a, b and x respectively. The negativity is 
invariant under interchange of b and x. For clarification see figure [?!). To calculate the 



X 



Figure 7. The chain separated into three parts: the impurity-spin (a), the part whose 
entanglement with the impurity-spin is to be quantified (b) and the part that separates 
the two (x). 

negativity between the impurity-spin and a part of the chain that is separated from the 
impurity-spin by a different part, we have to compute the reduced density- matrix that 
results from tracing out the separating part of the chain. We will now explain how this 
is done within our variational framework. 

Let the result of the generalized eigenvalue problem be the variational ground-state 
iPg given in the variational basis: ipg = 'Y^a^ipi. To take trace of the density-matrix over 
the region x in figure [TJ we express the density-matrix in a product-basis formed by 
states of region x and states of the rest of the chain. This requires us to represent the 
variational ground-state ipg in this product-basis. In doing so we follow the approach 
of reference [lOj, where the entanglement-entropy was calculated from a variational 
ground-state. For a chain with an odd number of spins these bases will in the following 
be introduced. We start with the case that the region x contains an even number of 
sites followed by the case that the region x contains an odd number of sites. For a chain 
with an even total number of sites an analogous calculation was performed. 
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6.1. An even number of sites to be traced out 

(i) For the impurity-spin we choose the 5'2-basis and denote the two states by /i = ^ 

and f2=i ■ 
(ii) For the region x we denote the basis by Qi. We take the region to contain at least 

two sites. The states are chosen to be: 

p n— 1 dimcrs 



91= XI °" t - ■■■ - t 



"■- ^ 2R sites 

r, n—1 dimcrs 



92= ^an t - ■■■ - i 



"- ^ 2R sites 

93 = t - ■ ■ ■ -t , 94 = i - ■ ■ ■ -i , 

2R sites 2R. sites 

^5 = ; - • • • -t , 96 = t - ■ ■ ■ -i , 

2R sites 2R sites 

2_R sites 

where the symbol '— ' is used for valence-bonds between neighbouring sites and R 
is the highest possible number of dimers that can could be formed in region x. 

(iii) For region b we define the states in the following way: 

L-2R~1 sites 
iY„i n—R—1 dimers 



h2 = 


2^ an t - • • • - !-•••- 

"=■^+1 L-2R-1 sites 
^_-|^ n—R—1 dimers 


h3 = 


2_^ an i - ... - !-•••- 



"—^+1 L-2R-1 sites 

In order to express the density-matrix of the impurity-spin and the region b in 
an orthonormal basis it is necessary to orthonormalize the basis {hi} as it is not 
orthogonal the way it is defined above: 

{h\h3)= E ^n{-2)-^^^U^] . (29) 

n=R+l Vv2y 

We do this following the usual Gram-Schmidt procedure. Our orthonormalized 
version is given by: 

h 



\/Chi\hi] 



h 



2 



\J{h2\h2 
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hs-{h3\hi)hi 
h-i = 



13 

(30) 



'{h\hs)-\{hs\h^)\^ 

For the product-space of the impurity-spin and region b we use the basis-states: 
h = fihi, k2 = fih2, h = fihs, 

h = f2hl, h = /2/l2, h = f2h3 ■ 

With these definitions, the variational ground-state can be written as: 



(31) 



where the matrix C is given by: 
/ 



c 








V 



V2 






llfelll 














2 











2 


aol / 


\\h2\\ 
V2 














K 

2 


























\\h2\\ 
2 


















/ 



where k = y H/^aP — |(/i3|/ii)P and \\h\\ = ^^/{h\h). 

With the matrix Gij := {gi\gj), whose matrix-elements can be computed using 
the prescription given in section [5l one can write the reduced density-matrix of the 
impurity-spin and region b as 

[Pabjij ^ v^^x Paxbjij 

^ / _^^in^jm\9n\Qm) i W'^j 



or in matrix form: 



Pab — CGC 



6.2. An odd number of sites to he traced out 



(33) 



(i) For the impurity-spin we choose the S^-basis and denoted the two states by /i = '|" 

and f2=\r ■ 
(ii) For the region to be traced out we denote the basis by gi. The states are chosen to 

be: 

9i= - ■■■ -i , 92 = - ■■■ -i , 



2R+1 sites 
R 



2R+1 sites 



n— 1 dimers 



^3 = X^ fln t - 



t - • • • -t , 



n=l 



2R+1 sites 



Impurity Entanglement in the J — J2 — 6 Quantum Spin Chain 



14 



R 



n—1 dimcrs 



9i 



95 



y^Qn i - 



- t- 



-t> 



n=l 
R 



2R+1 sites 
n—1 dimcrs 



E«- L_ 



t 



4., 



71=1 

R 



2R+1 sites 
n—1 dimers 




t - • • • -{ , 

V 

2R+1 sites 
5 5'8 = j— ■ ■ ■ — , 
2_R+1 sites 2_R+1 sites 

where the symbol '—'is used for valence-bonds between neighbouring sites and R 
is the highest possible number of dimers that can could be formed in region x. 
(iii) For region b we define the states in the following way: 

/ii = t- • • • - , h2 = i- 



L-2R-2 sites 



L-2R-2 sites 



h. 



N-l 



n—R—1 dimers 



t 



"=■^+1 L-2R-2 sites 

We orthonormalize and build a product-basis for the impurity-spin and region b 
just as it is done in the case of an even number of sites that have to be traced out. 

With Gij := {gi\gj) and 
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where k = y (h^lh^) — |(/i3|/ii)P, we can again write the reduced density-matrix as 

Pab = CGC^ . (34) 

For a chain with an even total number of sites we followed an analogous procedure. In 
order to save space we have not included it here. 

6.3. The Negativity 

With the expression for the reduced density-matrix of the impurity-spin and region b, 
Pab, at hand it is now straight forward to calculate the negativity according to its 
definition [26] 

EIA.I-1 



Af 



(35) 
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6.4- Precision of the Variational Approach 

The variational states that form our starting point should yield almost exact results 
close to the MG-point ( J2 = J/ 2) where a chain with an even number of spins is fully 
dimerized. However, the states we include in our variational calculation (see section [3]) 
form only a subset of the states in which the bulk of the chain is dimerized. At the 
MG-point the contribution arising from these neglected states is negligible and the 
variational approach is very good. However, as J2 is decreased toward the critical point 
J| we expect the contribution from the neglected states to grow in importance and at 
J|, where the ground-state is no longer dimerized, we expect our variational approach 
to fail. Nevertheless, for intermediate J2, J| < J2 ^ J/2 we still expect the variational 
approach to yield qualitatively good results. 




>.0.3 



40 60 

site 



Figure 8. Variational results for the negativity (filled symbols) together with DMRG- 
data from reference [llj (empty symbols) for two values of the impurity-coupling Jj^. 
The calculations were done for a chain of length L — 250, with frustration J2 = 0.42 
and dimerization 6 = 0. 

It is possible to calculate the negativity. A/", using DMRG techniques, as was shown 
by Bayat et al [llj. In figure [8] we compare our variational results for A/" to DMRG 
results from reference ^T] at J2 = 0.42 J for 2 different values of Jk'- Jk = 0.16 and 
Jk = 0.4 (Here (5 = 0). The definition of A/" used in reference [11] differs from the one 
used by us by a factor of two. To account for this difference we multiplied their data by 
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two. For both values of Jx we observe a very good agreement between the DMRG and 
variational results. 

7. J\f Compared to Simp at the MG-Point 

Having established the validity of the variational approach we now turn to a comparison 
of the two measures for the impurity entanglement A/" [11] and Si^p [9\ . In figure |9] we 
show results for both A/" and Si^^p at the MG-point for a system with L = 200 sites 
and an impurity-coupling Jk = 0. Shown is also the impurity entanglement arising 
from a single impurity valence bond (IVB) [9j between the impurity-spin and the bulk 
of the chain. For Jk = and L even, the idea of an impurity valence bond is easy to 
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Figure 9. Variational results for the negativity JV and the impurity entanglement 
Sirnp scaled by l/(21n(2)) at impurity-coupling Jk = and 6 = 0. Shown is also 
the impurity entanglement arising from an impurity valence bond (IVB) scaled by 
l/(21n(2)). The calculations were done for a chain of length L = 200 at the MG-point. 
For clarity only a fraction of all data points are shown. 



understand. The uncoupled impurity must form a singlet with the unpaired S* = 1/2 in 
the bulk of the chain. If the entire system (impurity and bulk) is divided in 2 parts A 
and B this impurity valence bond will contribute either ln(2) or depending on whether 
the unpaired spin is in part B or A. With the wave-function for the unpaired spin (psoi{x) 
and P = Jj^ \fso\{x)\'^dx, the impurity entanglement is then simply [9l[T0]: 



S; 



IVB 
imp 



;i-p)ln(2). 



(36) 
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To a first approximation ipso\{x) can be taken to be the wave-function of a free particle 
in a box [10]: 



y^so\{i) 



m 
sm , 
L \L 



(37) 



It is also possible to obtain more precise estimates of yjsoi [10] • Apart from an overall 
scaling factor of 2 ln(2) the two measures are in the case shown in figure M essentially 
identical. 
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(b) L = 199 



Figure 10. Variational results for the negativity JV and the impurity entanglement 
Sirnp scaled by l/(21n(2)) at impurity-coupling Jx = 1- Tlie calculations were done 
at the MG-point for a chain of length L ~ 200 (a) as well as L = 199 (b). For clarity 
only a fraction of all data points are shown. 



We now compare the two measures of the impurity entanglement in the limit 
Jk = 1- In figure [TOTa) we show variational results for both quantities at Jk = 1- Since 
the impurity here is coupled with the same strength as the remaining sites and since 
we consider a chain with an even number of sites, the chain is in its dimerized ground- 
state. The dimerization of the state implies that every spin is maximally entangled with 
one of its neighbours. It follows that the impurity is only entangled with its neighbour 
and not with any other spin. While this can be inferred from the negativity shown in 
figure [TOT a) which is exactly for / > 0, this is not the case for S'imp which in fact 
becomes negative (!). The difference in the two measures arises from the subtractive 
procedure used to define S'imp where the term S'[/(l,/ — 1,L — 1) in addition to terms 
cancelling out Su{l,l,L) also contains a contribution arising from the single unpaired 
spin which is present because L — 1 is odd. This contribution from Su{l, I — 1,L — 1) 
can be identified as the single particle entanglement [10] (SPE). Hence, Simp fails in this 
case to be a good measure of the impurity entanglement since S[/(l, / — 1, L — 1) is not 
a good reference state and cannot be identified with S'(no impurity). It would be very 
interesting to define a measure of the impurity entanglement based on a relative entropy 
but again one would need a well defined reference state with the impurity absent. 
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For a system with an odd number of sites the SPE is the only contribution to 
'S'imp at Jk = 1. Here the impurity is viewed in the general sense as arising from the 
boundary conditions yielding the non-zero S'imp- However, the negativity, since it is 
really a tri-partite measure, is not sensitive to this and is negligible for / > 1 as shown 
in figure [TOTb) at the MG-point. 

It may be conjectured that JV and Simp agree well as long as the IVB contribution 
dominates Simp as is the case in figure ID We also note that the SPE vanishes at Jg ilQ] 
and we expect A/" and S'imp to exhibit the same scaling behaviour at J| as it indeed has 
been demonstrated [11 [11]. For J2 > Jg, A^ is the more faithful measure of the impurity 
entanglement and in the following we mainly use this measure. 



8. The Negativity and the SPE for general Jk at the MG-Point 
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Figure 11. (a) Variational results for the negativity for different Jk with L = 200. 
(b) Variational results for the negativity for fixed size, I = 40, of the subsystem x as a 
function of Jk with L = 200. All results are obtained at the MG-point with 6 = 0. 



As we saw in the previous section, the impurity entanglement for an even length 
system with J^ = can be seen as arising from an impurity valence bond (IVB). 
At the MG-point for Jk = this was calculated within the variational approach in 
reference [10]. Here we focus on the negativity A/" and the complete cross-over as Jk is 
varied between and 1. Variational results of A/" at different Jk for a system of length 
L = 200 (including the impurity site) are shown in figure [TTT a). At J^- = the impurity 
entanglement is long-range, extending throughout the chain, but as Jk is turned on it 
quickly becomes suppressed and for Jk larger than ~ 0.3 it has all but disappeared 
reflecting the fact that the impurity valence bond now preferentially terminates close to 
the impurity site. For these non-zero values of Jk the appropriate (fsoi to use would then 
no longer be the wave-function of a free particle, as in equation (137|) . but instead a wave- 
function describing a localized state bound to the impurity site. The presence of such a 
localized state should be reflected in an exponentially decaying impurity entanglement 



Impurity Entanglement in the J — J2 — 6 Quantum Spin Chain 



19 



away from the impurity site. For Jk > 0.15 we have verified that the negativity TV 
shown in figure [TTT a) indeed has a tail decaying exponentially with the site index. The 
length scale, ^loc, associated with this exponential decay is therefore clearly distinct from 
the bulk spin-spin correlation length in the system which here is effectively zero. The 
sharpness of the cross-over is illustrated in figure [TTT b) where for a chain of 200 spins 
the negativity with the region x containing 40 spins is plotted versus J^- As can be 
seen, JV transitions from ~ 0.5 to very close to Jk = 0.15. For Jk < 0.15 we have 
not been able to identify any exponentially decaying part in A/" as obtained within the 
variational approach even for systems substantially longer than 200 and it appears that 
a finite J^ ~ 0.15 is needed to induce an exponential decay in JV. This would imply 
that the presence of a non-zero Jk cannot be modelled as a simple one-dimensional 
potential well, since this would always have a bound-state independent of the depth of 
the potential well. Likely, an enlarged variational space for the calculation will change 
the value for J^ and could possibly drive it all the way to zero. However, Bayat et 
al |JJJ, using DMRG, do not find an exponential decay for J2 = 0.46 and Jk = 0.16 
consistent with a non-zero J^. Due to the complexity of the calculation of A/" in larger 
sub-spaces we have, however, been unable to complete such calculations. 
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Figure 12. Variational results for the S'imp for different Jk with L = 199. All results 
are obtained at the MG-point with 6 — 0. For clarity only a fraction of all data points 
is shown. 



In a similar manner we can study the single particle entanglement (SPE) for general 
Jk within the variational approach. Our results are shown in figure [121 The single 
particle entanglement arises from the presence of a single unpaired spin for odd L 
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which is always present in the ground-state for odd L. The unpaired spin moves freely 
throughout the system when Jk = 1 but as Jk is decreased away from 1 towards 
zero it rapidly becomes localized on the impurity site thereby quenching the impurity 
entanglement. As can be seen in figure [T2] the impurity entanglement rapidly transitions 
around Jk ~ 0.8 — 0.85. While the variational approach agrees closely with DMRG data 
away from this transition region, some quantitative differences appear in this region due 
to the limitation of the Hilbert space used in the variational approach. The variational 
and DMRG results are however qualitatively the same and we have for clarity not 
included DMRG results in figure [121 We have verified that for J^ < 0.85 the data 
for Simp shown in figure [12] develop an exponentially decaying part with an associated 
length scale, ^loc, different from the bulk spin-spin correlation length. Also here, it 
appears that ^loc diverges before Jk reaches 1 but this effect could possibly be due to 
finite-size effects in the DMRG calculations or limitations in the sub-space used in the 
variational calculations. 

9. Impurity Entanglement with 5^0 
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Figure 13. (a) Variational results for ipsoi (markers) compared to Ai(i/^ + z) (lines) 
(see main text), (b) Variational results for the negativity at impurity-coupling Jk — 
for three values of the dimerization S plotted along with the scaled IVB entanglement 
(lines). The calculations were done for a chain of length L — 500, with J2 = J/2. 



The variational approach that we employ here is straight forward to use for any 
J2, 6 and Jk and we now discuss our results for the impurity entanglement for non- 
zero explicit dimerization 6^0. As we stressed in the introduction the impurity 
entanglement can be long-range in these systems due to the near degeneracy of the 
two singlet ground-states of the chain with periodic boundary conditions. With a non- 
zero 6 this degeneracy is lifted and we expect the impurity entanglement to decrease 
dramatically. Essentially, this is because when 6 is increased it is more costly for the 
impurity to be entangled with far away parts of the systems: If the impurity is in a 
singlet with a spin far away, there are more dimers on weak bonds than if the impurity 
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bonds with a site that is close to it. If we consider L even, then, in order to have 
sizeable entanglement, the bulk of the chain must have a site in a valence-bond with the 
impurity-spin and this spin is then bound to one end of the chain pE] . 

For a chain with an even number of spins we expect to be able to describe TV in 
terms of an IVB picture for which an estimate of ipsoi is needed for non-zero 6. Such 
an estimate has been derived by Uhrig et al (reference [IS]). They showed that to good 
approximation ipso\{i) oc Ki{i/^ + z), where Ai(x) is the Airy-function and zi = —2.3381 
is its biggest root. Here, the size of the area the impurity-spin binds to is roughly given 
by ^ = (3r7i5/2)^^/^, where m ^ l/(l + 7/-\/65). Our variational approach directly yields 
the lowest state vector which is simply (/^soi- In figure [TST a) we compare this directly 
to the Ai{i/C, + z) where it can be seen that the agreement is good. It is noteworthy 
that the maximum in (fsoi is quite distant from the end of the chain in agreement with 
previous results. 

We can use this result to estimate A^ using equation (136|) . In figure [T3r b) we 
see typical graphs of the behaviour of the negativity as 6 is increased with J^ = 
and L = 500. The range over which the impurity-site is entangled is decreased upon 
increasing the dimerization. Figure [T^ b) also includes the impurity valence bond 
entropy (IVB) obtained from equation (l36l) using the Airy functions for (/^soi- For easier 
comparison with the negativity we scaled by a factor of 1/2 In 2. For 5 = we used 
equation (1371) for ipsoi- A reasonable agreement between A^ and the IVB contribution 
(solid lines) is observed. Overall, as one would expect, the impurity entanglement is 
now dramatically suppressed. Even for S = 15 x 10~^ the impurity entanglement is 
effectively zero beyond 60 sites. Following the above discussion we expect the extent 
over which the impurity entanglement is non-zero to diverge as 5^^^^ as dictated by ^. 

10. Conclusion 

We have presented detailed discussion of how variational calculations can be used to 
calculate the negativity. This approach is quite generally applicable and should be very 
precise close to the MG-point. A complete characterization of the impurity entanglement 
for any Jk can then be obtained. For J^ = {L even) the impurity entanglement is 
long-range and extends throughout the system while for larger Jk > 0.15 it decays 
exponentially with a length scale clearly different from the correlation length in the 
system. For L odd the SPE present at J;^ = 1 is rapidly destroyed once Jx is decreased 
below Jk = 0.8 — 0.85. In both cases it would be valuable to perform higher precision 
calculations in order to determine whether or not a critical Jk or (1 — J^) is needed 
to destroy the impurity entanglement. When a non-zero dimerization is introduced and 
the ground-state degeneracy is lifted, the impurity entanglement rapidly disappears. It 
would be interesting to investigate impurity entanglement in other gapped systems with 
degenerate singlet ground-states. 



Impurity Entanglement in the J — J2 — S Quantum Spin Chain 22 

Acknowledgments 

We are grateful to I. Affleck and N. Laflorencie for interesting discussions and prior 
collaborations that initiated this work. We also acknowledge fruitful discussions with 
A. Bayat, S. Bose and P. Sodano and thank them for making their DMRG data for the 
negativity available to us. This research was supported by NSERC and CFI and was 
made possible by the facilities of the Shared Hierarchical Academic Research Computing 
Network (SHARCNET:www.sharcnet.ca). 



References 

[1 
[2 



[3 
[4 
[5 
[6 

[7 

[8; 

[9 
[10 

[11 
[12, 
[13 
[14 

[15 
[16 

[17: 



Calabrese P and Cardy J. Entanglement entropy and quantum field theory. J. Stat. Mech., 

P06002, 2004. 
White S R. Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett., 

69:2863, 1992. 
Zhou H-Q, Barthel T, Fjaerestad J, and SchoUwock U. Entanglement and boundary critical 

phenomena. Phys. Rev. A, 84:050305(R), 2006. 
Laflorencie N, S0rensen E S, Chang M-S, and Affleck I. Boundary effects in the critical scaling of 

entanglement entropy in ID systems. Phys. Rev. Lett., 96:100603-4, 2006. 
Cho S Y and McKenzie R H. Quantum entanglement in the two-impurity Kondo model. Phys. 

Rev. A, 73:012109, 2006. 
Kopp A and Le Hur K. Universal and measurable entanglement entropy in the spin-boson model. 

Phys. Rev. Lett., 98:220401, 2007. 
Le Hur K, Doucet-Beaupre P, and Hofstetter W. Entanglement and criticality in quantum impurity 

systems. Phys. Rev. Lett, 99:126801, 2007. 
Le Hur K. Entanglement entropy, decoherence, and quantum phase transitions of a dissipative 

two-level system. Annals of Physics, 323:2208-2240, 2007. 
S0rensen E S, Chang M-S, Laflorencie N, and Affleck I. Impurity entanglement entropy and the 

Kondo screening cloud. J. Stat. Mech., LOlOOl, 2007. 
S0rensen E S, Chang M-S, Laflorencie N, and Affleck I. Quantum impurity entanglement. J. Stat. 

Mech., 2007:P08003, 2007. 
Bayat A, Sodano P, and Bose S. Negativity as the entanglement measure to probe the Kondo 

regime in the spin-chain Kondo model. Phys. Rev. B, 81:064429, 2010. 
Sodano P, Bayat A, and Bose S. Kondo cloud mediated long-range entanglement after local quench 

in a spin chain. Phys. Rev. B, 81:100412, 2010. 
Eriksson E and Johannesson H. Impurity entanglement entropy in Kondo systems from conformal 

field theory. Arxiv preprint arXiv:1102.2492, 2011. 
Affleck I, Laflorencie N, and S0rensen E S. Entanglement entropy in quantum impurity systems 

and systems with boundaries. Journal of Physics A: Mathematical and Theoretical, 42:504009, 

2009. 
Bose S. Quantum communication through an unmodulated spin chain. Phys. Rev. Lett., 

91:207901, 2003. 
Christandl M, Datta N, Ekert A, and Landahl A J. Perfect state transfer in quantum spin networks. 

Phys. Rev. Lett., 92:187902, 2004. 
Christandl M, Datta N, Dorlas T C, Ekert A, Kay A, and Landahl A J. Perfect transfer of 



Impurity Entanglement in the J — J2 — S Quantum Spin Chain 23 



[18 

[19 

[20 

[21 

[22 

[23 

[24 
[25 
[26 
[27 

[28 

[29 

[30 

[31 

[32 

[33 

[34 

[35 

[36 

[37 

[38 

[39 

[40 

[41 

[42 

[43 



arbitrary states in quantum spin networks. Phys. Rev. A (Atomic, Molecular, and Optical 

Physics), 71:032312, 2005. 
Burgarth D and Bose S. Perfect quantum state transfer with randomly coupled quantum chains. 

New J. Phys., 7:135, 2005. 
Burgarth D, Giovannetti V, and Bose S. Efficient and perfect state transfer in quantum chains. 

J. Phys. A, 38:6793, 2005. 
Wojcik A, Luczak T, Kurzynski P, Grudka A, Gdala T, and Bednarska M. Unmodulated spin 

chains as universal quantum wires. Phys. Rev. A, 72:034303, 2005. 
Plenio M B and Semiao F L. High efficiency transfer of quantum information and multiparticle 

entanglement generation in translation-invariant quantum chains. New J. Phys., 7:73, 2005. 
Campos Venuti L, Degli Esposti Boschi C, and Roncaglia M. Long-distance entanglement in spin 

systems. Phys. Rev. Lett, 96:247206, 2006. 
Campos Venuti L, Degli Esposti Boschi C, and Roncaglia M. Qubit teleportation and transfer 

across antiferromagnetic spin chains. Phys. Rev. Lett., 99:060401, 2007. 
von Neumann J. Gott. Nachr., 273, 1927. 

Wehrl A. General properties of entropy. Rev. Mod. Phys., 50:221, 1978. 

Vidal G and Werner R F. Computable measure of entanglement. Phys. Rev. A, 65:032314, 2002. 
Affleck I, Kennedy T, Lieb E H, and Tasaki H. Rigorous results on valence-bond ground states in 

antiferromagnets. Phys. Rev. Lett, 59:799-802, 1987. 
Affleck I, Kennedy T, Lieb E H, and Tasaki H. Valence bond ground states in isotropic quantum 

antiferromagnets. Commun. Math. Phys., 115:477-528, 1988. 
Tribedi A and Bose L Quantum critical point and entanglement in a matrix-product ground state. 

Phys. Rev. A, 75:042304, 2007. 
Tribedi A and Bose L Entanglement and fidelity signatures of quantum phase transitions in spin 

liquid models. Phys. Rev. A, 77:032307, 2008. 
Hirano T and Hatsugai Y. Entanglement entropy of one-dimensional gapped spin chains. J. Phys. 

Soc. Japan, 76:074603, 2007. 
Path G, Legeva O, Lajko P, and Igloi F. Logarithmic delocalization of end spins in the s = f 

antiferromagnetic Heisenberg chain. Phys. Rev. B, 73:214447, 2006. 
Fan H, Korepin V, and Roychowdhury V. Entanglement in a valence-bond solid state. Phys. Rev. 

Lett, 93:227203, 2004. 
Alipour S, Karimipour V, and Mermarzadeh L. Entanglement and quantum phase transitions in 

matrix-product spin-1 chains. Phys. Rev. A, 75:052322, 2007. 
Geraedts S D and S0rensen E S. Exact results for the bipartite entanglement entropy of the AKLT 

spin-1 chain. Journal of Physics A: Mathematical and Theoretical, 43:185304, 2010. 
Majumdar C K. Antiferromagnetic model with known ground state. Journal of Physics C: Solid 

State Physics, 3:911-915, 1970. 
Haldane F D M. Spontaneous dimerization in the S = 1/2 Heisenberg antiferromagnetic chain 

with competing interactions. Phys. Rev. B, 25:4925, 1982. 
Eggert S. Numerical evidence for multiplicative logarithmic corrections from marginal operators. 

Phys. Rev. B, 54:9612, 1996. 
Eggert S and Affleck L Magnetic impurities in half-integer-spin Heisenberg antiferromagnetic 

chains. Phys. Rev. B, 46:10866, 1992. 
Laflorencie N, S0rensen E S, and Affleck L The kondo effect in spin chains. J. Stat Mech., P02007, 

2008. 
Shastry B S and Sutherland B. Excitation spectrum of a dimerized Next-Neighbor 

antiferromagnetic chain. Phys. Rev. Lett., 47:964, 1981. 
S0rensen E S, Affleck I, Augier D, and Poilblanc D. Soliton approach to spin-Peierls 

antiferromagnets: Large-scale numerical results. Phys. Rev. B, 58:R14701, 1998. 
Gaspers W J and Magnus W. Some exact excited states in a linear antiferromagnetic spin system. 

Phys. Lett A, 88A:103, 1982. 



Impurity Entanglement in the J — J2 — S Quantum Spin Chain 24 

[44] Gaspers W J, Emmett K M, and Magnus W. The Majumdar-Ghosh chain. Twofold ground state 

and elementary excitations. J. Phys. A, 17:2687, 1984. 
[45] Beach K S D and Sandvik A W. Some formal results for the valence bond basis. Nuclear Physics 

B, 750:142-178, 2006. 
]46] Uhrig G S, Schonfeld F, Laukamp M, and Dagotto E. Unified quantum mechanical picture for 

confined spinous in dimcrized and frustrated spin chains. The European Physical Journal B, 

7:67, 1999. 
]47] Sutherland B. Systems with resonating-valence-bond ground states: Gorrelations and excitations. 

Phys. Rev. B, 37:3786, 1988. 



